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PREMIÈRE  THÈSE. 


EXPOSITION 

DE  LA 

MÉTHODE  DE  M.  GYLDEN 


POUR  LE 

DÉVELOPPEMENT  DES  PERTURBATIONS  DES  COMÈTES. 


INTRODUCTION. 

1.  Dans  un  Mémoire  couronné  par  l’Académie  des  Sciences  de  Paris 
( Mémoire  sur  le  calcul  des  perturbations  qu  éprouvent  les  comètes , Supplé- 
ment aux  Comptes  rendus,  t.  I),  Hansen  a proposé  en  1847,  pour  le  dé- 
veloppement analytique  des  perturbations  des  comètes,  une  méthode 
nouvelle,  consistant  essentiellement  à partager  l’orbite  de  la  comète 
en  plusieurs  parties,  et  à former  des  développements  particuliers  pour 
chacune  d’elles.  Malgré  l’accroissement  de  convergence  qui  en  résulte, 
la  méthode  de  Hansen  n’a  pas  été  appliquée,  la  multiplicité  des  déve- 
loppements diminuant  singulièrement  les  avantages  que  la  méthode 
semblait  avoir,  et  l’on  continua  à calculer  les  perturbations  d’année  en 
année  par  la  méthode  des  quadratures.  En  juin  1869,  M.  Hugo  Gylden 
publia,  dans  le  Bulletin  de  l' Académie  des  Sciences  de  Saint-Pétersbourg, 
un  Mémoire  sur  le  même  sujet.  M.  Gylden  reprend  l’idée  de  Hansen,  et 
réussit  par  un  artifice  ingénieux,  en  introduisant  les  fonctions  ellip- 
tiques, à augmenter  la  convergence  à la  fois  par  rapport  aux  deux  va- 
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riables  dont  dépendent  les  positions  de  la  comète  et  de  la  planète  per- 
turbatrice. Ce  premier  Mémoire  de  M.  Gylden  fut  suivi  de  plusieurs 
autres  insérés  dans  les  Recueils  des  Académies  de  Saint-Pétersbourg, 
de  Copenhague,  de  Stockholm,  et  relatifs  à la  théorie  des  fonctions 
elliptiques.  Ces  derniers  Mémoires  renferment  un  grand  nombre  de 
méthodes  et  de  résultats  remarquables,  et  intéressent  autant  les  géo- 
mètres que  les  astronomes. 

2.  Nous  nous  sommes  proposé,  dans  ce  travail,  d’exposer  aussi  sim- 
plement et  aussi  brièvement  qu’il  nous  a été  possible  de  le  faire  la 
nouvelle  méthode  de  calcul  des  perturbations.  Les  idées  de  Hansen 
étant  peu  répandues,  nous  en  avons  fait  une  exposition  générale.  Nous 
avons  terminé  ce  Mémoire  par  un  exemple  numérique  destiné  à mon- 
trer dans  quelles  conditions  la  méthode  de  M.  Gylden  peut  s’appliquer 
avantageusement.  Pour  y parvenir  plus  clairement  nous  avons  choisi 
l’exemple  qui  avait  été  traité  par  Hansen  dans  le  Mémoire  cité  plus 
haut,  c’est-à-dire  le  calcul  des  perturbations  que  la  Terre  exerce  sur 
la  comète  d’Encke,  et  nous  avons  divisé  l’orbite  de  la  comète 
comme  l’avait  fait  Hansen  lui-même.  On  verra  que,  dans  ces  condi- 
tions, la  méthode  de  M.  Gylden  ne  s’applique  pas  : l’orbite  doit  être 
divisée  en  un  plus  grand  nombre  de  parties,  de  telle  façon  que,  dans 
chacune  d’elles,  la  distance  des  deux  astres  varie  surtout  en  raison  du 
déplacement  de  la  planète.  Néanmoins  il  n’est  pas  nécessaire  de  multi- 
plier beaucoup  le  nombre  des  points  de  division  pour  que  la  méthode 
puisse  être  appliquée;  cinq  ou  six  pourront  suffire,  et,  comme  on  aug- 
mente beaucoup  la  convergence  des  développements  par  rapport  aux 
deux  variables , l’avantage  de  la  méthode  ne  peut  être  contesté. 

3.  Soient 

x,  y , z les  coordonnées  rectangulaires  d’une  comète  par  rapport  à des 

axes  passant  par  le  centre  du  Soleil; 
x',y',  z'  celles  d’une  planète  perturbatrice; 
m'  la  masse  de  cette  planète; 

A la  distance  de  la  planète  à la  comète; 

/■'  celle  de  la  planète  au  Soleil. 
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Les  équations  du  mouvement  de  la  comète  sont 
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[}.  désignant  un  nombre  constant  pour  tout  le  système. solaire,  et  il  la 
quantité 

, / i xx' rr'+  zz'\ 

,n  u — ^ — ); 

on  a supposé  nulle  la  masse  de  la  comète. 

A cause  de  la  petitesse  de  la  quantité  fi,  on  a coutume  d’intégrer 
d’abord  les  équations  précédentes  en  y supposant  O nul;  puis  on 
regarde  les  intégrales  du  mouvement  elliptique  ainsi  obtenu  comme 
représentant  encore  le  mouvement  vrai,  ce  qui  est  toujours  possible  si 
l’on  regarde  les  constantes  introduites  par  l’intégration  comme  des  fonc- 
tions du  temps.  On  peut  même  faire  en  sorte  que  non-seulement  x , y, 
z,  mais  aussi  leurs  dérivées  premières,  conservent  les  mêmes  formes 
dans  le  mouvement  vrai  que  dans  le  mouvement  elliptique.  On  obtient 
ainsi  six  équations  entre  les  dérivées  par  rapport  au  temps  des  élé- 
ments du  mouvement  elliptique,  ces  éléments  eux-mêmes  et  le  temps. 

4-,  L’intégration  des  équations  du  mouvement  elliptique  conduit  à 
introduire,  au  lieu  des  coordonnées  rectangulaires  x , j,  z,  les  coor- 
données polaires  dans  le  plan  de  l’orbite,  et  à remplacer  la  variable 
angulaire  dans  ce  plan  par  l’anomalie  moyenne.  Pour  les  besoins  de 
l’Astronomie,  on  est  ainsi  amené  a conclure  des  variations  des  éléments 
du  mouvement  elliptique  celles  du  rayon  vecteur,  de  l’anomalie 
moyenne  et  de  la  latitude  ou  plutôt  de  son  sinus;  mais  il  est  clair  qu’il 
n’y  a là  que  des  transformations  faciles,  et  que  la  véritable  difficulté 
de  la  méthode  réside  dans  la  formation  des  dérivées  de  fi  par  rapport 
aux  éléments  du  mouvement  elliptique. 

Le  premier  terme  de  la  fonction  perturbatrice  étant  beaucoup 
plus  grand  que  le  second  quand  la  distance  des  deux  astres  est  très- 


(J 


B.  BAILLAUD. 


petite,  et  ses  dérivées  partielles  dépendant  de  c’est  la  formation  de 

cette  quantité  qui  doit  présenter  le  plus  de  difficultés,  et  c’est  à son 
étude  que  nous  nous  bornerons  dans  ce  travail. 

A dépendant  des  anomalies  moyennes  de  la  comète  et  de  la  planète, 
ou  de  variables  équivalentes,  se  développera  en  une  série  périodique 
de  sinus  et  de  cosinus  de  multiples  de  ces  variables.  Comme  A2  est  de 
la  forme 

r2+r'2  — Arr'cosj/cos/'—  Bit'  sin/sin/'  — C rr'  cos/sin/'  — D rr’  sin/cos/', 

où  r et  r'  désignent  les  rayons  vecteurs  de  la  comète  et  de  la  planète, 
/et/'  leurs  anomalies  vraies,  on  doit  d’abord  former  les  quantités 

r-,  rcos/,  rsin/, 
r’\  r'cos/',  /•' sin/'. 

Il  est  évident  d’ailleurs  que  la  connaissance  de  ces  six  quantités 
équivaut  à celle  des  positions  de  la  comète  et  de  la  planète  dans  leurs 
orbites,  et  par  conséquent  doit  suffire  pour  former  toutes  les  quantités 
dont  dépend  le  calcul  des  perturbations. 

I. 

5.  Soient 

a le  demi-grand  axe; 
e l’excentricité  de  l’orbite  de  la  comète; 
u son  anomalie  excentrique; 
t le  temps. 

On  a les  relations  bien  connues 

/ /•  = a(i  — e cos «)» 

6 r cos/=  a(cos«  — e), 

1 | /■  sin/^-  a y i — t^sin//, 

f ndt  = - du , 

a 

i a ( i — e2  ) 

t i-t-ecos/ 

2 1 / J 

I r 2 

f n dt  — df. 

ci 2 y'  i — a" 
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6.  Soient  rt  et  r2  les  rayons  vecteurs  de  deux  points  de  l’orbite  situés 
de  part  et  d’autre  du  grand  axe.  11  est  aisé  de  représenter,  à l’exclusion 
de  tous  les  autres,  les  points  de  l’ellipse  compris  entre  eux  dans  la  por- 
tion qui  contient  le  périhélie.  Il  suffit  de  former  une  fonction  pério- 
dique d’une  variable  auxiliaire  qui  ait  pour  maxima  r,  et  r2  et  pour 
minimum  a( i — e).  Posons 

/’=  A + B sin#  -+-  C sin2#. 


Cette  fonction  aura  un  minimum  a[i  — e)  si  l’on  peut  l’écrire 

(3)  r — a[ i — e)  -f-  (M  sin#  H-  N)2, 

la  valeur  absolue  de  M n’étant  pas  inférieure  à celle  de  N. 

r sera  maximum  quand  cos x sera  nul,  c’est-à-dire  quand  x aura 
l’une  des  valeurs  — 90°  et  -+-  90°.  Ces  valeurs  maxima  seront  r2  et  r,  si 

“M  — N = \jr2  — a ( 1 — e), 

M -f-  N = \Ji\  — ci[i  — e) . 

Les  deux  radicaux  ayant  le  même  signe,  puisque  M2  est  plus  grand 
que  N2,  on  peut  les  prendre  positifs. 

La  comparaison  de  la  relation  (3)  avec  la  première  des  relations  (1), 
mise  sous  la  forme 

, , . 11 

r = a i — e M-  2 ae  sin2  - 5 
2 

donne 

,,,  . u f M . N \,2 

( 4 ) sin2  - = - — : sin.r  -+-  j • 

2 \V2  ae  \liaej 


Le  calcul  de  M et  N sera  facilité  si  l’on  pose 


ce  qui  donne 

tang(45°  — 9) 


M 

y/2  ae 


N 

= e cos<p,  — — = esinç, 
y 2 ae 


— 2 a ( 1 — e 
4 ae 


Si  l’on  prend  s positivement,  cos<p  sera  positif,  et,  par  suite,  9 sera 
compris  entre  — 90°  et  -+-  90°;  comme  tang(45°  — 9)  est  positive,  9 est 
compris  entre  — 45°  et  -f-  45°. 
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Les  relations  précédentes  prennent  des  formes  utiles  quand  on  y in- 
troduit les  anomalies  excentriques  a,,  u2,  qui  correspondent  aux  points 
/•,,  r.2.  On  a 

, \ / . "7 

yr2  — a[i  — e ) = 2.11e  sin  — 5 

\Jr , — ali  — e)  = iaes\n  — ■> 


puisque  et  l~  sont  moindres  que  180  degrés. 


On  en  conclut 


u , . Il, 

e cos©  = il  sin h sin  — ] > 

2 2 


, , M|  . U7 

s sin  9=7  sin .sin  — 


et,  par  suite, 


1 / . M, 

sin2©  = — sin2  — 
2 £2  \ 2 


sin2—  = 


>'  1 — i'i 

4 atz1  ' 


. U7 

sin  — 

iang(45°  — 9)  = , 

. **1 

sin  — 

2 


(5) 


H sin2  — -4-  sin2 


11  est  aisé  d’exprimer  r,  rcosf,  rsinf,  ndt  en  fonction  de  x en  par- 
tant de  la  relation  (4).  Nous  remarquerons  d’abord  qu’en  extrayant  les 
racines  carrées  des  deux  membres  de  cette  relation  on  doit  écrire 


(6) 


sin  - —dz  e(cos9  sin^-i-  sin<p), 


en  adoptant  le  signe  supérieur  quand  u est  moindre  que  ut  et  le  signe 
inférieur  quand  u est  plus  grand  que  u2,  car  la  quantité  entre  paren- 
thèses s’annule  au  périhélie,  et  en  s’annulant  change  de  signe.  Elle  est 
d’ailleurs  positive  pour  x — 90°,  car 

, , . U, 

e (cos  9 -h  sin  9)  = sin  — 
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D’après  cela,  les  formules  (i)  donnent  immédiatement 

r = a[ i — e -h  eeJsin’cp)  -+-2ae£2sin2tp  sin.r  — aee1  cos2©  C0S2Æ', 
r cos  f — a(  i — e — £2sin2cp)  — 2«E2sin2ç  sin^e  -f-  r<£1oos2ç  co.S2.r, 

rsin  f = 2 ne  \J i — e‘[  eoscp  sin  x -+-  sin  9)  \j i — £2sin2cp  — £2sin  2 9 sin x — e2cos29  sin ‘x, 


2 £ COS  © COS  X d.T 


ndt  = - 


\/ 1 — î2sin2cp  — £2  sin  2 cp  sin  x — £2cos2cp  sin2jr 

En  tenant  compte  de  l’observation  relative  à la  formule  (6),  on  recon- 
naît aisément  que  le  radical  qui  entre  dans  les  deux  dernières  formules 
est  positif. 

7.  L’expression  de  r en  fonction  de  l’anomalie  excentrique  étant 
analogue  a celle  de  ^ en  fonction  de  l’anomalie  vraie,  on  entrevoit  aisé- 
ment la  possibilité  de  représenter  l’arc  d’ellipse  qui  joint  les  points  rt 
et  r2  en  passant  par  l’aphélie.  Il  faudra  poser 


- = A'+  B'  sin.r  -t-  C'  sin2.r  ; 


r 


r aura  pour  maximum  a(t  -+-  e ),  si  cette  valeur  de  - peut  s’écrire 


-t-  ( M' sin.r  -+-  N')2, 


r a(i+e) 


en  supposant  que  la  valeur  absolue  de  M'  ne  soit  pas  moindre  que  celle 
de  N',  et  pour  minima  r,  et  r2 , si 


les  radicaux  étant  encore  pris  positivement. 


En  identifiant  la  valeur  précédente  de  - avec  celle  que  donne  la  pre- 


r 


mière  des  formules  (a),  on  trouve 


(4  bis) 
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On  facilitera  ici  le  calcul  de  M'  et  de  N'  en  posant 
M' 


■V- 


-*2) 


2 e 


= £'côsç',  N 


V- 


- — e si n o , 


ce  qui  donne 


e — r. 


lang(45"  — 9’)=  \/ ~ ~ 

, /,  , n[r,-h  r2)(n-  e)  — n\v. 

£ li-e) 


4 eri  r. 


les  radicaux  étant  pris  positivement,  de  sorte  que  ç>'  est  compris  entre 
— 45  et  -f-  45°. 

Par  l’introduction  des  anomalies  vraies /,  et f2  des  points  r , et  r2,  les 
formules  précédentes  se  simplifient.  On  a 


/ r i / 2 

a(n-é>)~~  V fl(i  _ e»)  cos— ’ 


/ I I / 2«  f 

V r>  a(i  + e)  y «(i  — e2)  C°5  2 ’ 

les  signes  des  seconds  membres  étant  déterminés  par  la  considération 
f f 

de  ceux  de  cos  — et  de  cos  — • Par  suite, 

2 2 


, , i / f fi\ 

e coscp  = - ( cos  ~ — cos  )• 

, . , i ( f , 

£ sinœ  — cos  - — h cos  — 

2 \ 2 2 


et  encore 


sm  2 ?'  — f cos2  ^ — cos2  rr  ) = rzz,rî~  «(>  — 


2 £ 


2 ) l\ez'‘‘rir-1 


f 

cos  — 

tang(45°—  ?')  = 

cos  — 
2 


/. 


[5  bis)  £,2=  ^ |^cos2  ^ + cos2 

Pour  déterminer  ^ cosy,  sin/et  ndt  en  fonction  de  x,  on  remarque 
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que  la  relation  (4  bis)  donne 

f 

cos  ^ =±  e'  ( cos cp'  sin.r  — sincp')  ; 


i i 


mais  on  doit  choisir  le  signe  supérieur,  car  si  la  parenthèse  change 
de  signe  quand  on  passe  à l’aphélie,  il  en  est  de  même  du  premier 
membre,  et,  pour  x = 90°,  la  parenthèse  est  positive. 

On  a donc 

f 

cos  - — £'(cos<p'  sin.r  — sin<p')> 

et  l’on  en  conclut 


a(  j — e2)  - = 1 — e -4-  ez'2-V-  c£'2sin2cp'  — 2/?£,2sin2ç/sin.r  — CE^eos’ç' C0S2.r, 

cos / = 2 -+-  £,2sin29'  — 1 — 2£'2sin2cp'sin.r  — e'2  cos2©'  cos2.r, 

2£'(cos©'sin^r—  sin  ©')  y'  1 — £'2sin2©'  -4-  e'2 si n 2cp'sin  x — £'2cos2cp'  sin2.r, 
r2  2 e'  coscp'  cosxdx 

a-\J  1 — e2  y/ 1 — s'2sin2cp'  -f-  £,2sin2<p'sin.r  — E^cos'cp'  sin2.*/ 

le  radical  qui  entre  dans  les  deux  dernières  formules  étant  positif. 

8.  Les  formules  précédentes  se  simplifient  beaucoup  quand  l’un  des 
points  de  division  est  situé  au  périhélie  ou  k l’aphélie,  et  les  formules 
ainsi  obtenues  permettent  de  représenter  l’ellipse  partagée  en  quatre 
portions.  Par  une  transformation  facile,  on  peut  faire  en  sorte  que  ces 
quatre  parties  soient  obtenues  en  faisant  varier  la  variable  de  zéro  à 
36o  degrés.  Rangeons  ces  parties  dans  l’ordre  suivant  lequel  on  les 
rencontre  en  partant  du  périhélie  et  tournant  dans  le  sens  du  mouve- 
ment de  la  comète.  Les  formules  relatives  à la  première  et  à la  quatrième 
se  déduiront  de  celles  du  n°  6,  dans  lesquelles  x varie  de  — 90°  à 4-  90°, 
ou  si  l’on  veut  de  4-  270°  k 4-  90°,  puisqu’il  n’entre  pas  dans  ces  for- 
mules de  sous-multiple  de  x;  pour  la  deuxième  et  la  troisième  partie, 
on  appliquera  les  formules  du  n°  7,  où  x varie  de  + 90°  k 4-  270°. 
On  fera  les  substitutions  suivantes  : 


sin  f = 
ndt  — 


Première  et  deuxième  région  : x = 2 y — 90°, 
Troisième  et  quatrième  région  : x = 2 y -f-  90% 
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et  l’on  obtiendra  successivement  les  quatre  régionsen  faisant  varier  y de 
zéro  à 90  degrés,  de  90  à 180  degrés,  de  i3o  à 270  degrés,  de  270  à 
060  degrés. 

Nous  désignerons  par  E , , E2<  E3,  E,  les  valeurs  de  z ou  de  z dans  ces 
quatre  régions,  par  4>(,  <î>2 , <I>3,  O.,  celles  de  y et  de  9'.  Les  formules  sont 
alors  les  suivantes  : 


Première  région  : 

. , . . Jr — ai  1 — e)  . u,  .11  _ . 

= Ei=  \ =:  sin  , sin  - = E,  sm*r, 

y 2 ae  2 2 

r — a{\  — e + -f  eEJ)  — aeYL\  cos  2/  -f-jaeE1 2 *  cos4/, 

r cos  f = a[i  — e — f E2)  + a E?  cos 2/  — ~a EJ  cos 4 /, 

/■sin/  = aE,  \/i  — e*  (1  — cos  2/)  \J  1 — E2  sin'/, 

, r 2E1  sin  2r 
n dl—  - J 

Deuxième  région  ; 

4V=-4  5°,  E,=  i/ 


a sj  1 — EJ  sin4/ 


dy. 


. a(i  -h  e)  — r,  f f . 

1 — e)  — - = cos  — ■>  cos  - = E2sm2r, 

22 


1er, 


a ( 1 — e 
r 


1 — e -h  |eEj  — eE2  cos2/  -4- -j-eE2  cos 4/, 

cos/=:fEJ2  — 1 — E,  co s 2/  +|E2cos4/, 
sin  f — E,(i  — cos  2/)  \J  1 — E2  sin4/, 

ndl  = — 


/■2  2E,  sin2r  , 

---  — - — dy. 

«2y/i  — e2  yi  — E2sin4/ 


Troisième  région  : 


<!>.  =+  45°,  E3  = ( 1 — e) 


-4-  e ) — r2 


2e/’; 


f f 

cos  — j cos  - =—  E3sin’K, 
22 


n ( 1 — e2 


= 1 — e + f E2  — eE3  COS2/  -f-  jeE2  cos4/, 

cos/=  jE2  — 1 — E2cos2/4-jE2cos4/, 
sin/=  E3(cos2/ — 1)  ^1  — E2sin4/, 


dt  = 


r2  2E3sin2/ 

a}\J\  — e2)  1 — E2 sin4/ 


dy. 
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Quatrième  région  : 
$,=  -45»,  E* 


V 


r2  — fl  i — r)  . u j . u 

— = sin  — 5 sin  - = E4sin2r, 

2 lie  2 2 


r = a[  \ — e 4-{eEJ)  — ne  EJ  cos2/  4-  \ae EJ  cos4/> 
rcosf  = a( i — e —{EJ)  4-  «EJcos2/  — {«EJ  cos4/. 


r si n/  = «E,  ^ i — e2  (cos 2/  — i ) 1 — EJ  sin4/, 
r 2E4  sin2r 


ndt  = 


a \j  \ — E4  sin1/ 


dy. 


9.  Dans  les  valeurs  de  r,  rcos/,  r sin/,  ndt  données  dans  les  nos  6 
et  7,  les  termes  périodiques  sont  multipliés  par  s2  et  s'2.  Il  est  manifeste 
que,  si  ces  nombres  s2  et  s'2  sont  très-petits,  les  développements  des 

quantités  rsin/,  ndt , en  séries  procédant  suivant  les  multiples 

de  x,  convergeront  rapidement.  D’autre  part,  les  valeurs  de  s 2 et  s'2 
données  par  les  formules  (5)  et  (5  bis ) montrent  clairement  que  les 
quantités  E2,  E2,...  sont  moindres  que  les  quantités  £2,  s'2,  de  sorte 
que  les  formules  du  n°  8 devront  être  employées  quand  celles  des  nos  6 
et  7 ne  conduiront  pas  à des  développements  assez  convergents. 


10.  Il  est  intéressant  d’examiner  comment  varient  les  quantités  E2, 
E2  quand  le  point  r{  occupe  sur  l’ellipse  toutes  les  positions  possibles. 
On  reconnaît  de  suite  qu’à  mesure  que  ce  point  s’éloigne  du  périhélie 
E,  augmente  et  E2  diminue.  Il  est  donc  impossible  de  faire  décroître 
ces  deux  quantités  en  même  temps;  il  peut  être  utile  de  rechercher  le 
minimum  de  la  somme  E2  4-  E2. 

Si  l’on  pose 

e = sin41, 

on  a les  formules  bien  connues 


>4 


d’où  l’on  lire 


/ 
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tang  - = tang  7 + -1-  long  -> 


cos2  - 


4 2 


7T  dit 


COS2  - 


tang  , , , - - 

4 2 

2 2 

rdf=  cos  <4  du. 


L’expression  + E^,  égale  à 
a pour  dérivée 


. u f 

sin2 — 1-  cos’  -■> 
2 2 


. U U . f f 

sin  - cos  - du  — sin  - cos  - df. 

2 2 2 2 

En  remplaçant  sin  ^ cos  " et  sin  '^cos  par  des  quantités  proportion- 

nelles, et  de  même  du  et  df,  on  voit  que  le  signe  de  cette  dérivée  est 
celui  de 

r2 

COS2  O/. 

a-  T 

Il  s’ensuit  que  la  somme  considérée  est  minimum  pour 

r=  a cos  41. 

sin/  = sin»#,  f ~~ r — u, 


Alors 

et,  par  suite, 


cos-  = sin-  et  E2  = E2 
2 2 


De  plus  la  troisième  relation  citée  donne 


tang2  — = tang  + — ) > 


d’où 


sin2  — 1-  cos 
2 2 


tang  ~ ' 


On  voit  que,  pour  le  minimum,  chacune  des  quantités  E?  et  Ej  est 
moindre  que  et  que  ces  quantités  sont  d’autant  plus  petites  que  Tel- 
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lipse  est  plus  excentrique.  Les  mêmes  conclusions  s’appliquent  visible- 
ment aux  quantités  Eg  et  E2. 

11.  Ce  résultat  conduit  à rechercher  pour  quelles  positions  du 
point  r,  les  quantités  E2  et  E2  sont  moindres  que 

Il  faut  d’abord  que  u soit  un  angle  aigu,  à cause  de  la  valeur  de  E, . 
La  relation 

cos2  — <i~ 

2 


est  satisfaite  toutes  les  fois  que  l’angle  f est  obtus. 

Si,  au  foyer  de  l’ellipse  d’où  parlent  les  rayons  vecteurs,  on  mène 
une  perpendiculaire  au  grand  axe,  elle  rencontre  l’ellipse  en  un  point 
dont  la  distance  au  foyer  est  a{i  — e2).  L’arc  d’ellipse  compris  entre  ce 
point  et  le  sommet  du  petit  axe  est  le  lieu  des  points  considérés. 


12.  Dans  le  cas  où  l’on  emploie  les  formules  des  nos  6 et  7,  il  est 
intéressant  de  rechercher  encore  quels  doivent  être  les  points  de  sépa- 
ration pour  que  la  somme  £2  + s'2  soit  minimum.  Comme 


£2  + 


r,  -t-  r,  — 2 a ( 1 — e) 
4 ae 


r,  et  r2  étant  indépendants  l’un  de  l’autre,  la  somme  considérée  sera 
minimum  quand  on  aura 

r,  = r2  = a \/ 1 — e1—  a cos^. 


Ce  résultat  ne  différant  de  celui  du  n°  10  que  par  la  situation  du 
point  r2,  la  plupart  des  calculs  indiqués  à cette  occasion  s’appliquent 
ici;  on  trouve 

. , , . , «2  , fl 

sin2  — = cos2—  = sm2—  = cos2  — 5 
2222 

de  sorte  que 

e2  — e'2  et  e2  4-  s'2  = 2 sin2—1  = 1 — lang  — • 

2 2 

Les  points  de  division  qui  correspondent  au  minimum  de  la  somme 
£2  -+-  s'2  sont  donc  les  mêmes  que  dans  le  cas  du  partage  de  l’ellipse  en. 
quatre  parties,  et  les  valeurs  minima  de  c2  -h  e'2  et  de  E2  + E^  sont 
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égales.  Il  ne  faudrait  pas  en  conclure  que  l’on  ne  gagne  rien  à employer 
les  formules  du  n°  8.  Celles-ci  ne  diffèrent  de  celles  des  nos  6 et  7 que 
par  un  changement  d’argument  qui  ne  modifie  vraiment  pas  les  for- 
mules et  par  la  substitution 

E2  = 2 £2, 

de  sorte  que,  dans  le  cas  du  minimum,  le  module  E2  des  formules  du 
n°  8 est  réellement  la  moitié  de  celui  des  formules  des  nÜS  6 et  7. 

13.  1!  peut  être  nécessaire  d’augmenter  beaucoup,  dans  une  portion 
peu  étendue  de  l’orbite,  la  convergence  des  développements.  On  y par- 
vient aisément  en  ne  représentant  que  cette  partie  de  l’ellipse.  Soient 
r3  et  /•,  les  rayons  vecteurs  de  deux  points  situés  d’un  même  côté  de 
l’axe.  Supposons,  pour  fixer  les  idées,  que  le  point  rt  soit  dans  la  pre- 
mière moitié  de  l’ellipse,  et  le  point  r3  entre  le  point  i\  et  le  périhélie. 

Les  formules  du  n°  8 permettent  de  représenter  l’arc  compris  entre 
le  périhélie  et  le  point  r{  au  moyen  d’une  variable  y variant  de  zéro 

à é-  Si  l’on  veut  ne  représenter  que  l’arc  r3r{,  il  suffit  d’exclure  les 

points  voisins  du  périhélie,  et,  pour  représenter  l’arc  compris  entre  le 
périhélie  et  r ,,  on  exclura  les  points  compris  entrer,  et  r3.  On  est  ainsi 
conduit  à poser 

sinj  = A sin  z,  cos/  = Icosii, 

avec  la  relation 

h ‘ -h  l2  — i , 

h et  l étant  les  sinus  et  cosinus  de  l’anomalie  y relative  au  point  /;3. 

En  faisant  varier  z de  zéro  à on  représentera  l’arc  compris  entre 

le  périhélie  et  r3,  et,  en  faisant  varier  u de  zéro  à on  représentera 

l’arc  r3r,.  Les  termes  périodiques  de  r,  rcos/,  ...  seront  multipliés 
par  des  puissances  de  / ou  de  h ; si  l’arc  r3r,  est  suffisamment  petit, 
/ sera  très-petit,  de  sorte  qu’on  aura  notablement  augmenté  la  con- 
vergence. 

On  trouvera  dans  le  Mémoire  cité  de  Hansen  les  formules  relatives 
à ces  anomalies  a,  que  Hansen  nomme  anomalies  intermédiaires.  Nous 
remarquerons  seulement  qu’elles  donnent  le  moyen  de  partager  l’el- 
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lipse  en  autant  de  parties  qu’on  voudra  et  d’augmenter  en  quelque 
sorte  indéfiniment  la  convergence  des  développements  en  séries. 


14.  La  théorie  exposée  dans  les  paragraphes  précédents  est  entière- 
ment conforme  aux  idées  de  Hansen.  L’introduction  des  fonctions  ellip- 
tiques permet  d’augmenter  encore  la  convergence  des  développements. 

Supposons  qu’on  ail  partagé  l’orbite  de  la  comète  en  deux  parties, 
et  que  les  points  de  division  soient  symétriques  par  rapport  au  grand 
axe.  L’anomalie  partielle  relative  a la  partie  inférieure  est  alors  définie 
par  l’équation 

.u  , 

sin  - = nr  s sin#, 

2 

où 


. « , 

£ = sin  — ? 

2 

et  l’on  a 

r — a ( i — e ) 4-  2 ae  e2  sin2  #, 
r cos f = a(i  — e)  — 2«s2  sin2#, 


r sin  f — iaz\J  i — e2  sin#  \Ji  — £2  sin2#, 

, (i — e 4- 2e£2sin2#)  2£  cos#  , 

ndt  — — dx, 

s!  i — £2sin2# 


et  les  formules  relatives  à la  partie  supérieure  sont 


où 

et 


cos  — — d sin#, 
2 


= i — e 4-  2ê£'!  sin  2#, 

r 

cos f = — i 4-  2e'2  sin2#, 

sin/  = 2e'  sin#  ^ i — e'2  sin2#, 

, r2  2e'cos# 

ndt  — — - — - - ■ 

a2\/ 1 — e2  \J i — e'2  sin2# 


dx. 


3 
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15.  Nous  poserons  pour  la  partie  inférieure 

2 E 

x — am — co  (mod.e), 


E désignant  l’intégrale  elliptique  complète  de  première  espèce  relative 
au  module  a.  Alors  les  valeurs  de  r,  rcos/,  ...  deviennent 


2 E 

r = a( i — e ) H-  2«e£2sin2am  — - co, 

TC 

„ . , . 2 E 

r cos/  — a[\  — e)  — nav  sin2  am  — - co, 

TT 

2.E  . 2 E 


rsin/*=  2 cï s i/i  — e2sin am  — co  Aam  — co, 

TT  7T 

2 E 


mc//  = — ( i — e + 2e£J  sin2  am  — - m ) 2£  cosam  — co.rfco, 


2 E 


et  ces  valeurs  peuvent  être  développées  en  séries  très-convergentes 
procédant  suivant  les  sinus  et  cosinus  des  multiples  de  co. 

On  a d’abord 

2 E . 2 E 7T  (1  2 E 

sinam — co  Aam  ■ — co  = — cosam  — co. 

TC  TC  2 E d CO  TC 

Le  facteur  ndt  exige  le  développement  de 

. , 2E  2 E 

stn’am- — -co  cosam  — co. 

TC  TC 


Posons 


2E  TC  . . 

am  — - w = - — amj  ( mod.  e)  ; 


l’expression  proposée  deviendra 

cos2amj  sinam/. 

C’est  une  fonction  uniforme  doublement  périodique;  elle  a les  mêmes 
périodes  que  sinamj;  elle  a deux  infinis  simples  et  deux  infinis 
doubles;  elle  change  de  signe  quand  y augmente  de  E.  C’est  donc  une 
fonction  linéaire  de  sinamj  et  de  sa  seconde  dérivée,  c’est-à-dire  de 
2 e 

cosam  — co  et  de  sa  dérivée  seconde.  En  calculant  cette  seconde  déri- 
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2 E 2 E 1 / 7r  \ 2 c/2  2E  2 E 

e’ sin2am  — m cosam — m = - — ^ -7—  cosam- h | cosam  — m. 

TC  TC  2 \ 2 E / «M1  71  7T 


Les  formules  relatives  a la  partie  inférieure  deviennent  ainsi 

, . , . , 2 E 

r = a[  1 — e)  + 2ae£2  sin2am  — co, 

7T 

rcosf  — a ( 1 — e)  — 2«£Jsin2am  — m, 

7T 


rsin/  = — 2 rt  \j  1—  e1  ( — ^ ^ cosam  — co, 


7T  \ (1 

2 E / c/cO 


2E 

TC 


dt  /2E\  2E  t c \ d2  2E 

ri  — 2 — s cosam  — co  h-  2e  — = z -7—  cosam  — co. 

d CO  V 7T  / TC  \ 2 L / C/CO2  7T 


• • • 2 Îj  2 Ij 

On  sait  d’ailleurs  développer  sin2am  — co,  cosam  — co,  et  par  suite  sa 

dérivée  seconde  en  séries  toujours  convergentes  procédant  suivant  les 
cosinus  des  multiples  de  co.  On  a effectivement 


2E  2 E . , — / 1 a „ (p 

£ cosam  — co=  — 4 v^/  cos  co  h - — cos  3 m -t-  — - — r cos  5 m 

TC  TC  1 \ X -t-  Q 


1 -+-  qi 

2E  0 / TT  \ 2 / q 2<72  3qJ 

e’sin’am — m = 8 — H — - -+-  - — 

r.  \2E/  \i  — q2  i—q*  1 — qb 


1 -+-  qi 


8(^ 


2 E ) \ 1 — q 


q 2 e2  . 3 q3  c 

1 C0S2M-1 — COS4co  + --■  -COSOM- 


ou 


E, 


q = e 


E,  étant  l’intégrale  elliptique  complète  de  première  espèce,  relative  au 
module  s,  complémentaire  de  e. 


16.  Le  développement  des  formules  relatives  a la  partie  supérieure 
exige  plus  d’attention. 

Posons  encore 


x — am 


(mod,  £*), 


3. 


E'  étant  l’intégrale  complète  de  première  espèce  relative  au  module 


et  a)  variant  de  - à 3 -• 
2 2 

Nous  aurons 


a i — el 


, ■ 2 E' 

i — e 4-  2ee 2 sm2am  - — - « 


cos  f — — i -+-  2 s'2  sin2am 


2E' 


. „ n d 2E' 

sin  / = — 2£  — — -7-  cosam  • — w, 

47  O IV  TT 


idt  — 


2E'  c/m 

2 e'  2 E'  r2  2 E'  , 

■ ■ ■ — • - — ■ — cosam wdo, 


de  sorte  que  la  véritable  difficulté  réside  dans  le  développement  de 
qui  peut  s’écrire 


ali+e 


2e 


n e'2sin2am w 


2 E' 


i — e 


7T 


On  peut  procéder  de  la  manière  suivante. 

Soient  

e',  = \A  — e'2 

et 


alors 


y/— = tangamv  (mod.  e',); 


a{l~'re)  i -+-  e'2tang2amv  sin2am  o 


ce  qui  peut  s’écrire 


al  i + e 


= î 


2E' 

e'2  tang2amv  sin2am  ■ — - m 

7T 

TT  2Ë^ 

i + e'2tang2amvsin2am  w 


Remplaçons  tangam(v,  e',)  par 
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nous  aurons 


aE' 


e'ssin2am(v  i , z)  sin2am  — — m 


a ( i + e i 


i + 


2 E' 


i — e^sin’am  (v — i,  z' ) sin2am « 

TT 


ce  qui  peut  s’écrire 

/’  _ tangam  (v\J—  i , z'] 

a[i-\~e)  A am  [v\!~  i , z') 


aE' 


e'^sinam  [v\J—\,z')  cosarn  [v\/—  i,  e')  Aam  (v y/—  i,  z')  sin2am  — 


X 


i — s'2 sin2am  ( y \J—  i , z')  sin2am 


2 E' 


et  d’après  des  formules  connues 


,,  s , , , dW  v \j—  î , « 

r i -i-  ^ ~ sinam(v>  ei)  cosam(v>  ei)  V 71  ) 


2E' 


a{  i -t-  ej 


A am  ( v,  e',  ) 

, / 2 E'  \ 

''hH 

la  fonction  n ayant  la  signification  que  lui  attribue  Jacobi. 

Les  deux  facteurs  du  dernier  terme  peuvent  être  développés  en  séries 
procédant  suivant  les  sinus  et  cosinus  des  multiples  de  « : on  utilisera 
pour  le  premier  la  relation 


sinam(v,E,)  cosam  (v,  e, 
Aam(v,  e',)  • 


£3  I;  loe  A am  (»,*', 


17.  Le  développement  de  r dans  le  dernier  cas  peut  encore  s’effec- 
tuer par  la  formule  du  binôme.  En  remplaçant  dans  la  valeur  de  r 
2 E' 

sin2am  — « par  son  développement,  on  aura 

r i 

a(i-t-e)  A -f- B^r  C0S2&) -t- Cç2cos4«  + • . • ’ 

A,  B,  C étant  des  coefficients  numériques. 

Le  second  membre  de  cette  égalité  peut  s’écrire 


i 

A 


i 


COS  2 O)  -+- 


R 


« 
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OÙ 


Si  l’on  pose 


on  aura 


R = ~ (C^Jcos4&)  -t- . . .). 


i 


i-f-  — q cos 2 ai 
A 7 


a i -4-  e) 


i 


i 


B i + cp 

i + t a cos  2 w 
A 7 


A cause  de  la  petitesse  de  9,  le  dernier  facteur  pourra  être  développé 
en  une  série  très-convergente.  Le  précédent  et  son  carré  qui  entre  dans 
la  valeur  de  n dt  peuvent  l’être  par  différentes  méthodes.  Si  l’on  pose 


ce  facteur  devient 


B _ 2 a 

A ^ 1 -+-  a2’ 


1 + a 2 

1-4-2 a COS20)  -4-  a2’ 


dont  le  développement  est 


1 + a2 
1 — a 2 


(1  + 2 a cos  2 w -l-  2a2  cos4co  -4- . . 


I 


son  carré  est 


(i  + a5)a 

(1  — a2)3 


[1  -4-  «2 4-  2 et  cos 2 a>  -4-  (3  — o2)a2cos4w 


18.  Les  formules  développées  dans  les  trois  derniers  paragraphes 
permettent  de  développer  toutes  les  quantités  nécessaires  à la  détermi- 
nation de  la  position  de  la  comète  en  des  séries  très-convergentes  pro- 
cédant suivant  les  multiples  des  anomalies  partielles,  en  supposant  que 
l’orbite  ait  été  divisée  en  deux  parties,  les  points  de  division  étant 
symétriques  par  rapport  au  grand  axe. 

Comme  on  le  verra  dans  la  suite  de  ce  travail,  une  telle  division  de 
l’orbite  est  très-souvent  insuffisante.  Il  sera  très-facile  d’introduire  de 
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nouveaux  points  de  division  sans  recommencer  le  calcul.  Si  l’on  veut 
en  placer  un  au  périhélie,  il  suffira  de  poser 

I — COS2W, 

sin cô  — sin2aj,  : • 

2 


En  donnant  à w,  les  valeurs  de  zéro  à -■>  on  n’obtiendra  que  la  moitié 
de  la  partie  inférieure  qui  commence  au  périhélie.  Posons  de  même 

. , I — C0S2&), 

s i n m = — sin2  Gh  — •> 

2 

et  faisons  varier  w2  de  — ^ à zéro,  nous  aurons  représenté  l’autre  moitié 
de  la  partie  inférieure  de  l’orbite. 

La  partie  supérieure  de  l’orbite  ne  présentera  pas  plus  de  difficultés. 
Si  l’on  veut  placer  un  point  de  division  à l’aphélie,  on  fera  les  mêmes 

substitutions,  et  l’on  feravarier  oo,  de  ^ à n pour  la  première  moitié,  et 
de  n à ~ pour  la  seconde. 

S’il  est  encore  nécessaire  de  partager  les  portions  inférieures  de  l’or- 
bite, on  posera,  comme  il  a été  dit, 

COS  M|  = l COS  U, 

sinco,  = h sin  z, 
où 

A2 -4-  l7—l, 

et  l’application  répétée  de  telles  substitutions  permettra  de  diviser  l’or- 
bite en  autant  de  parties  que  l’on  voudra. 

II. 

19.  Dans  la  Section  précédente  on  a exposé,  d’après  les  idées  de 
Hansen  et  de  M.  Gylden,  le  moyen  de  représenter  par  des  séries  très- 
convergentes  procédant  suivant  les  sinus  et  cosinus  des  multiples  d’un 
certain  angle  les  quantités  r,  rcos/,  rsin /,  ndt.  11  est  aisé  de  tirer  de 
ces  développements  celui  de  A2  au  moyen  des  considérations  suivantes. 


2/j 
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Envisageant  chacune  des  planètes  et  des  comètes  du  système  solaire 
comme  décrivant  une  ellipse  suivant  les  lois  de  Kepler,  on  aperçoit  de 
suite  que  la  connaissance  de  leurs  positions  initiales  et  du  mouve- 
ment de  l’une  d’elles  conduit  à la  connaissance  du  mouvement  de 
toutes  les  autres. 

Désignons  par  co  l’anomalie  partielle  pour  la  portion  de  l’orbite  de  la 
comète  que  nous  voulons  étudier,  par  M l’anomalie  moyenne  de  la  co- 
mète, par  M'  celle  de  la  planète,  par  n'  son  moyen  mouvement  et  posons 

M'  = c + n'  t, 

c'  désignant  une  constante  dépendant  de  la  position  initiale  de  la  pla- 
nète et  changeant  de  valeur  à chaque  révolution  de  la  comète. 

Si  l’on  écrit 

, n' 

Al  — c —t-  — ni. 

n 

on  voit  qu’un  calcul  bien  simple  donnera  le  développement  de  M'  en 
une  série  procédant  suivant  les  sinus  et  cosinus  des  multiples  de  co. 
D’autre  part,  les  formules  de  mouvement  elliptique  permettent  d’en 
déduire  les  développements  de 

r',  r'  cos/',  r'sinf. 

Quelques  multiplications  permettront  enfin  d’en  déduire  le  développe- 
ment de  A2,  développement  qui  se  présentera  sous  la  forme  suivante  : 

A2  P o H—  Pi  si  n co  — i—  P2  c o s 2 co  —i—  P3  s i n 3 co  —h  . * . , 

P0,  P,,  P2,  ...  étant  des  séries  procédant  suivant  les  sinus  et  cosinus  de 
c'  et  de  ses  multiples. 

Nous  supposerons  que  les  points  de  divisions  de  l’orbite  aient  été 
tellement  choisis,  que  les  coefficients  des  polynômes  P,,  P2,  ...  soient 
très-petits  par  rapport  à ceux  de  P0  .11  faut  pour  cela  que  la  variation 
de  A2  dépende  surtout  de  la  variation  de  c'  et  très-peu  de  celle  de  co. 
L’exemple  qui  sera  développé  plus  loin  montrera  clairement  la  portée 
de  cette  condition. 

20.  On  a vu  que  le  calcul  des  perturbations  que  la  planète  produit 
dans  le  mouvement  elliptique  de  la  comète  dépend  surtout  de  la  for- 
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mation  de  la  quantité  Le  développement  en  série  de  cette  quantité 

est  la  partie  vraiment  difficile  du  calcul.  A2  est  représenté  par  une  série 
procédant  suivant  les  sinus  et  cosinus  des  multiples  de  c'  et  des  mul- 
tiples de  w.  Cette  série  a un  assez  grand  nombre  de  termes,  et  les  pre- 
miers sont  souvent  du  même  ordre  de  grandeur.  L’introduction  des 
anomalies  partielles  a permis  d’augmenter  la  convergence  relativement 
à la  variable  «.  Nous  allons  exposer  actuellement  le  procédé  ingénieux 
au  moyen  duquel  M.  Gylden  augmente  la  convergence  par  rapport  à c' . 
Soit  R l’ensemble  de  tous  les  termes  qui  suivent  P0; 

A2  = P0  4-  R, 

d’où 

a-’  = p7  — 4RP/  + ^R2P7.... 

Les  puissances  entières  successives  de  R s’obtiennent  aisément  par  la 
multiplication.  La  véritable  difficulté  réside  dans  la  formation  des  puis- 
sances de  P0. 

Soit 


P,  = A + A,  sine'  -h  A2sin2e'  + . . . + B,  cosc'  -4-  B2  cos2c'  + . . . . 


Les  coefficients  A2,  B2  seront  ordinairement  très-petits  par  rapport  à 
A,  A,,  B,.  Si  l’on  écrit 


P0  = A -+-  A,  sine'  + B,  cosc'  -t-  II, 


on  aura 

P0  — : (A  -+-  A,  sine'  -t-  B,  cosc') 


H V 

A -1-  A,  sine'  -+-  B,  cosc'/ 


Toute  la  difficulté  du  développement  des  puissances  de  P0  réside  dans 
celui  des  puissances  du  premier  facteur.  Posons 

A,  = — Msino,  B.^Mcosç,  — = y.. 

Le  premier  facteur  de  P0  devient 


A [1  -+-  cos  (c'  4-  9)]. 


4 


2.6 
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Si  la  dislance  des  deux  astres  peut  devenir  très-grande  dans  la  partie 
considérée  de  l’orbite,  le  coefficient  A sera  grand,  ce  qui  facilitera  la 
formation  des  puissances  négatives  de  P0.  Si  cetle  distance  peut  devenir 
petite,  [j.  sera  presque  égal  à l’unité,  et  le  développement  des  puissances 
négatives  de 

A [i  -+-  y.  cos  [c'  -+-  <p)] 

conduira  à une  série  très-peu  convergente. 

21.  Ce  facteur  peut  s’écrire 

A|,  + 'i)[,~7^sinHlc'  + F)]' 

Posons 

— h2,  \ (c'  4-  F)  = am  — \ (mod./r), 

1 -+-  [I.  7T 

K étant  l’intégrale  elliptique  complète  de  première  espèce  relative  au 
module  k.  Le  facteur  précédent  deviendra 

A (i  -+-  /a)  A2am  (mod./r). 

Ce  facteur  et  ses  puissances  pourront  être  développés  en  séries  très-con- 
vergentes de  sinus  et  de  cosinus  des  multiples  de  £. 

L’argument  £ devra  être  substitué  à l’argument  c' . La  substitution 
devra  être  effectuée  partout,  de  sorte  qu’il  sera  nécessaire  de  former 

les  développements  des  sinus  et  cosinus  des  multiples  de  am  — 2,  des 

2K 

puissances  des  sinus  et  cosinus  de  am  — Ces  questions  peuvent  être 

résolues  de  bien  des  manières.  Il  s’agit  ici  de  calculs  numériques  bien 
plus  que  de  développements  analytiques.  M.  Gylden  utilise  une  méthode 
appliquée  à d’autres  questions  par  Hansen  dans  plusieurs  de  ses  Mé- 
moires. 


22.  Soit,  par  exemple,  à développer 


Posons 
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K'  étant  l’intégrale  elliptique  complète  de  première  espèce  relative  au 
module  k'  complémentaire  de  k,  et 

• fn  2 K v(2n)  ^ v(2w)  > 

sin2nam — £ = X„ X,  coS2£ 

T.  0 i — q2  2 

-^xr>c°S4Ç-...-^x«cos2/ 

1 H I — Q 


D’après  la  relation  facile  à établir 


. . 2Kr 

a2sin2nam- — -£  . „ 

TT  ~ 2 K 

fTT : — = 211(211  — i ) sin2"-2  am  — £ 


— 4«J(i  + A2)  sin2nam  — £ -+-  2/1(2»+  1)  ^2sin2n+2am  — £, 

7T  ' TT 


on  aura  l’équation  récurrente 

4P 

\ zn  / 

-4-  4 ( 1 -+-  A-2)  Xj£2n)  — 2 « ( 2 » 1 ) k2  x^"+2)  ? 


TT  \ v(2n)  / \ v(2n — 2) 

— X.  =-2R  2/1-1  X. 

2 K / 2‘  v ' 21 


équation  qui  permet  de  calculer  de  proche  en  proche  les  coelïicients 
de  même  rang  dans  les  développements  des  puissances  successives  de 

2K  w 

sinam  — t. 

TT 

Posons  à cet  effet 

Y (2  n) 

— n 

Y(2«  — 2)  Z'2'*’ 

A2î 


l’équation  (a)  deviendra 


o = — 2/1  (2/1—  1)  +[^4/i2  ( i -4-  k1)  — 4P  J P™~  211  (2/i+  0 k7p,»pw 

Il  est  aisé  de  ramener  cette  équation  récurrente  à ne  plus  renfermer 

4. 
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qu’un  seul  coefficient  variable.  Posons  en  effet 

F2„  = 

l’équation  (a)  deviendra 


4«2(i  + a-2)  — 4«'2 


Posons  en  outre 
nous  aurons 


2 n -4- 1 

F2n  — p 1 n n2  * 2n  ~ PînPm+1- 


Pin  FsrtTÜLt/ij 


in  -+-  1 r ^ 

I : = CJj,,  — /j  1 2„  ï 2«-f-2  C7?nn5j,i-H2  , 


et  si,  pour  abréger,  on  pose 

■,  m -k-  1 . jp  p 

A2j|4-2  -T  2/1  ^ 2/<-t-2  > 

2 /l  — I 

l’équation  (a)  prendra  définitivement  la  forme 

l)  I = 7J5-2 n ~ ^2«-h2^2/iîS'2n+2  > 

et  sous  cette  forme  elle  permet  de  réduire  s>2„  en  fraction  continue 


TS5  jn  


^2  n- 


L’expression  de  F2«  donne 


F,  = — d’où 
1 -+-  l( 2 


'1  + A-2)2’ 


et,  par  suite,  s»  est  racine  de  l’équation  du  second  degré 


1 — z + 


h" 


La  solution 


i + A2)2 
. + A-5 


Z = O. 


k 


donnerait 
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2 K 

Les  coefficients  du  développement  de  sin2"am  — | augmenteraient  in- 
définiment avec  n.  On  doit  donc  adopter  l’autre  solution 

i -i- /i%  d’où  pn  — i. 

Adoptant  la  valeur  de  1+  k 2 pour  celle  de  ®2„+2,  correspondant  à 
une  valeur  très-grande  de  in  + 2,  on  remontera  de  proche  én  proche 
au  moyen  de  l’équation  (b),  et  l’on  obtiendra  ainsi  toutes  les  valeurs 
des  quantités  zs,  par  suite  celles  des  quantités p,  et,  comme  011  connaît 
X(9  et  on  en  conclura  toutes  les  quantités  X. 


23.  Si  la  valeur  1 + k 2 de  ne  paraît  pas  assez  approchée,  on 

peut  remarquer  que,  si  l’on  pose 


-/,„(«  + k ) —J  , 


fln  aura  pour  limite  l’unité.  D’autre  part,  l’équation  (b)  devient 
, . * r.  , 7..,  4'2  ( k y (2/t  + i)2  ^ 

( c ) * — fin  I 1 « / 1 \ Il  ) 1 ( \ \ f M+2  1 

L 4,t  \2*V/  2 //  ( 2 /l  -f  2 ) ^ 

d’où  il  suit  que f2n  diffère  de  l’unité  d’une  quantité  du  deuxième  ordre 
par  rapport  à l--  Posons 

fin  — 1 I -+"  J 2 ni 

et  négligeons  les  termes  dont  l’ordre  dépasse  le  second,  nous  aurons 

, , 4 7’  / 7t  V k 2 

f —h2 f — 4 \ n 

J in  n ./  *11  + 2 

Si  l’on  admet  que 


on  aura,  aux  termes  près  du  troisième  ordre, 

T *J+4*‘2( 


4/2 1 

i * \ 

4n’ 1 

yu) 

2h  + 2 

=/; 

/;«= 


7T 

2K 


2/1  1 y 1 — k2 


Posant  encore 


k'-h  4 i2 


fin  1 + 


2K 


2 « + i ; 


- k 2 


'fini 


3o 
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substituant  dans  l’équation  (c),  négligeant  les  termes  dont  l’ordre  dé- 
passe le  troisième,  et  supposant 

f in+i—J  in  > 

on  obtiendra  f2n  aux  termes  près  du  quatrième  ordre,  et  ainsi  de  suite. 


24.  Si  l’on  n’a  besoin  des  puissances  de  sinam  | que  pour  des 

valeurs  peu  considérables  de  n , il  suffira  de  modifier  l’emploi  de  l’é- 
quation récurrente  (a),  de  façon  à calculer  d’abord  les  coefficients  des 

2 K 

développements  des  premières  puissances  de  sinam  — 2,  pour  en  dé- 
duire ceux  des  puissances  plus  élevées.  On  posera 

X2  n — 2 

2 i 

—7l  — q?n-  J. 


L’équation  (a)  deviendra 


la')  o = — 2 n ( 2 n -t-  i ) h ' 


+ [4  «’ 


i -h  /rs)  — 4 é 


7T 

Ik 


— m{in 


Opérant  comme  au  n°  21,  on  posera 

„ 2n[2n  + i]/Ç 


4 n2  ( 1 + /r J)  — 4 


7T 

Ik 


puis 


P I 2/1  + I 

72/ï  vJ2n  P2rt>  P-2n  I 2 q xJ2/i— 2> 

/I  2 Tl  r“  J 


et  l’équation  récurrente  (a')  deviendra 

( b'  ) I = P2/1  p-7n—2  P m P'in—i‘ 

On  traitera  cette  équation  comme  celle  du  n°21,  avec  cette  seule 
différence  qu’on  en  tirera  p2„  en  fonction  de  /32„_2,  et  que  les  valeurs 
initiales  sont  connues  exactement 


d’où 


f*  0 = 0,  Po  = O, 
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Il  y a exception  dans  le  cas  où  i = o.  Alors 

i _ 14  g» _i_  G2 

7°  — ^TTi  f4 0 P <1  - — q — 3 j( 2 ^(2)  ' 

25.  Des  procédés  tout  à fait  identiques  aux  précédents  s’appliquent 

2 K 

au  développement  des  puissances  impaires  desinam-^-  S,  et  à ceux 

des  puissances  de  cosam  — 2 etde  Aam—  S. Nousnenousyarrêterons 
pas.  Il  nous  suffira  de  dire  que  l’on  part  des  formes  suivantes  : 

sin  3at  + . . . , 


sin2"-^1 

am 

2 K 

4_s/?_  x 

(2»+0 

si  n je  + 

4 vV 

X^0: 

TT 

1 - 7 

1 — q 3 

cos2" 

am 

2K 

1 

1 

T 

>* 

8q 

Y2n  cos 

2 AT  -f- 

i6r/2  . 

77 

® 1 

~ 7 

1-7’ 

cos2n+' 

am 

2k 

y 

(2«-H) 

COS  AT 

4vV 

yO-H-i; 

77 

i-q 

1 + q‘ 

) 3 

1 a m 

2 K 

? — 

: Z (*°  Hr  ■ 

87 

rj  (7«)  iÔfl2 

-Z  COS  2 AT  H ±— 

cos3at 


• (2-0 


tt  ' 0 i — q2  ‘ i — q“  • 

A^'am  — | = Z(2"+,)+  Z^,+,)  cos  2*  + Z(u+,)  cos 4 * +. . . 

tt  0 i -+-  q1  i 4-  q2  4 ^ 

et  des  formules 


^/2cosnamo7 

dx- 


— n ( n — 


k'2c  os"-1  ami 


-+-«’(  i — 2k'1)  cos"  ami  — n ( n -+-  1)  /î'2cos"+2  am  at, 


d 2 A"am  a: 
dx2 


n[n  — 1 ) /f'2  A"~2am  at  -t-  n2  (1  -4-  /r'2)  A"ami-;i(n  + i)  A',+2am  at. 


26.  La  même  méthode  peut  être  appliquée  au  développement  des 
fonctions 

2K  , 2K, 

sinraam- — 0 cosnam — £. 

TT  TT  ' 


Posons 


sin/i  am  — £ — 2 l\n)  sin ^ -{-  2l[n)  sin  2^  + 2 2(n)  sin  3?  -+- ... , 

7T  1 2 ^ 3 

cos/iam  ^ £ — r‘n)+  2r[n)  COS^  + 2r[";  COS 2;  + . . . . 
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On  a 


— (/2sinnam 


2k 


dp 

TC 

7 n2(i  4-  hH)  sin/i  am  £ 


2 K 


7 n ( n H-  i)  /r2sin  [n  4-  2)  am  — £ — 7 n ( n — 1)  /»2sin  (h 


2 am 


2K 


7T 

2 K 


1.  2k  . 

ircos/i  am  — £ 

TC 


dp 


2 K 


7 «2  (1  4-  h''1)  cos  n am  — £ 

TC 

{ « (n  4-  1)  /i2cos  (m  4-  2 ) am  2—  £ — 7 n [n  — 1)  /i2sin  ( n — 2)  am  — £ 

71  71 


On  en  conclut  les  relations  récurrentes 


| « ( n — 1 ) /r2  2 '"-2)  + ^n!(i4-  /r'2)  — i2 


in(u-i)  /.-* 4-  4 «2 f 1 4-  /, 


[* 


2k 

7 r 

2 K 


2‘“UiR(n4-i]K;"+,)  = o, 


4«-2)_1  I 1 ir|M)4-  i»(n  + I)*,I,(M  = 0, 


relations  que  l’on  pourra  traiter  absolument  comme  celles  des  numéros 
précédents. 

L’identité  de  forme  de  ces  deux  relations  montre  que 


de  sorte  que  le  rapport 


•^j.n-i-2)  <ç(n) 

pi»  T 2)  _ fW’ 
A i 1 

2(;n) 

ÏT3 


est  indépendant  de  n,  ou,  plus  exactement,  n’a  que  deux  valeurs  appli- 
cables, l’une  quand  n est  pair,  l’autre  quand  n est  impair. 

• 2 K 2 K 

Les  développements  connus  de  sinam  — § et  de  cosam  — \ et  de 
leur  carrés  montrent  que  les  deux  indices  d’un  quelconque  des  coeffi- 
cients ^ et  T sont  a la  fois  pairs  ou  impairs,  et  l’on  trouve,  en  faisant 
n — 1 ou  n = 2, 


y (2  n) 

^*9 


2,  ’ _ 1—  q11  2-li+i  '+  </2 

__  5 


1 1”}  I + T1'  ri::; 
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formules  qui  seront  de  la  plus  haute  importance  dans  le  calchl  des 
coefficients  2 et  r,  puisqu’elles  permettent  de  déduire  presque  immé- 
diatement les  coefficients  2 des  coefficients  r. 

Nous  terminerons  là  les  indications  relatives  au  développement  des 
sinus  et  cosinus  des  multiples  de  am  ^ £ et  des  puissances  des  fonc- 
tions elliptiques  élémentaires  de  £.  Dans  son  beau  Mémoire  ( Studien 
au f dem  Gebiete  der  Storungstheorie,  Mémoires  de  l’ Académie  impériale 
des  Sciences  de  Saint-Pétersbourg,  t.  XVI,  1871),  M.  Gylden  indique  de 
nombreuses  méthodes  pouvant  conduire  à la  formation  des  coefficients 
de  ces  développements.  Ses  travaux  sur  ce  sujet  intéressent  à la  fois 
les  géomètres  et  les  astronomes.  Nous  avons  essayé  de  réunir,  dans  ce 
travail,  ce  qui  nous  a paru  le  plus  immédiatement  utile  à ces  derniers. 
Nous  allons  actuellement  essayer  d’éclaircir  la  méthode,  en  en  com- 
mençant l’application  à un  cas  particulier. 

III. 

27.  L’objet  principal  de  cette  élude  étant  la  comparaison  de  la 
méthode  de  M.  Gylden  à celle  de  Hansen,  nous  avons  choisi  comme 
exemple  la  détermination  des  perturbations  que  la  Terre  exerce  sur  la 
comète  d’Encke,  ainsi  que  l’avait  fait  Hansen  lui-même  dans  son  Mé- 
moire couronné  ( Mémoire  sur  le  calcul  des  perturbations  qu  éprouvent 
les  comètes ).  Cherchant  surtout  dans  ce  calcul  un  élément  d’apprécia- 
tion, nous  ne  nous  sommes  pas  préoccupé  de  pousser  le  calcul  jusqu’au 
bout,  mais  seulement  de  former  le  développement  du  carré  de  la  dis- 
tance des  deux  astres.  Pour  cette  raison  aussi,  nous  avons  adopté  pour 
les  éléments  de  l’orbite  de  la  comète  et  de  celle  de  la  Terre  les  éléments 
dont  Hansen  s’était  servi,  savoir  : 

T = 1208",  1466, 
loga  = 0,3466760,  a' 

e = 0,8446760,  e' 

ra,  = i8ï«>48' 55', 8,  cj' 

Q = 334°29'  28",  8, 

/ = 1 3°  20'  40",  2 . 


= o , 0167796, 
= 1 oo1'  o'  1 3", o , 
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T désignant  la  durée  de  la  révolution  sidérale  de  la  comète,  Q la  lon- 
gitude de  son  nœud  ascendant  sur  l’écliptique,  sr,  l’arc  compris  entre 
le  nœud  ascendant  et  le  périhélie,  «l’inclinaison  de  l’orbite  sur  l’éclip- 
tique, a le  demi-grand  axe  de  l’orbite  de  la  Terre,  ë son  excentricité, 
rs'  la  longitude  de  son  périhélie.  Nous  avons  divisé  l’orbite  de  la  comète 
en  deux  parties,  les  points  de  séparation  étant  symétriques  par  rapport 
à son  grand  axe,  et  déterminés  par  cette  condition  que  leur  distance  au 
Soleil  soit  égale  à i.  Enfin  nous  n’avons  développé  les  calculs  que 
pour  la  partie  inférieure  de  l’orbite. 


28.  Désignant  par  a,  l’anomalie  excentrique  du  premier  point  de 
séparation,  par  u celle  d’un  point  quelconque  de  la  partie  inférieure 
de  l’orbite  et  par  co  l’anomalie  partielle,  on  a (n°  14) 


. u . i E . . . 

sm  - = rn  e sin — w mod.  e , 

2 TT  ' 

où 

. u , 

E = SUI  — , 

2 


co  devant  varier  de  — + le  signe  — se  rapportant  aux  valeurs 

de  co  comprises  entre  — ^ et  zéro,  le  signe  + à celles  qui  sont  com- 
prises entre  zéro  et  ‘j>- 
Posant  ensuite 

-E, 

g = -e~  h , 

E,  étant  l’intégrale  elliptique  complète  de  première  espèce  relative  au 
module  e,  complémentaire  de  e,  on  a trouvé 

log  £ = 1,621 1482, 

— = 24°  42'  24%  04, 

logE  = 0,2167170, 
logE,=  0,3652829, 
log  q — 2,07912t. 
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On  en  a conclu 

ecosam  — m — [1,61 58436]  cos  m 4-  [3,700144]  cos3<a 

TC 

-+-  [5,7793]  cos5m  -+-  [7,86]  cos 7 m -P  [8]  cosgoj, 

e2sinJam  — m — [2,g5i56i8]  — [2,941079]  cos2m 

— [3,321 168]  cos4m  — [5,5764]  cos 6m  — [7*78]  cos8m, 

et  ensuite,  par  l’application  des  formules  données  à la  fin  du  n°  14, 

r—  [ 1,8326789]  — [ i,5i5i47]  cos 2 m 

— [ 3,89524]  cos4m  — [4, i5o5]  cos6m  — [6,35 j cos8m, 

/•cos/ = — [2,718677]  -+-  [ i,588456]  cos 2 m 

4-  [3,968545]  cos4m  -4-  [4,2238]  cos6m  4-  [6,43]  cos8m, 

rsinf  — [ 1,9712033]  sinM  4-  [2,532625]  sin3M 

4-  [4,8336]  sin5M  4-  [5,o588]  sin7M  4-  [ 7,25]  singM, 

ni  — c = [ i,3o85i2i]  sinM  — 2,3i654i ] sin3M 

— [4,6622]  sin  5m  — [6,898]  sin  7&>  — [7,0]  singM  ; 

n est  exprimée  en  prenant  le  rayon  du  cercle  pour  unité  ; et  les  nombres 
entre  crochets  sont  les  logarithmes  des  coefficients  dont  ils  tiennent 
la  place. 

29.  Soient  zs  la  longitude  du  périhélie  de  la  comète,  f l’anomalie 
vraie  de  la  Terre.  On  a 

A:  = r‘  4-  r'1  — nrr'  sin2^cos[/'  4-/8-  [m'  — Q ) 4-  (ta  — Q)] 

- 2 rr'  cos3 ^ cos  [/'  / 4-  (ra'  - Q)  - (®  — Q)]. 

Les  nombres  donnés  au  n°  27  donnent 

(gj'  — Çl)  4-  (GJ  —&)==  3o8°  lÿ 

[m1  - Q)  - (ex--  Qj=  - 57018'  II", 6. 


5. 
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ou 


AJ  = r2  -4-  r'2  — A rr'  cos  f cos/'  — B/r'sin/sin/  ' 
— C/t'cos/ sin/'  — Drr'  sin/ cos /' 

logA  = 0,0344462, 
logB  = 0,0207989, 
logC  =0,2257077, 
logD  = 0,2146266, 


de  sorte  que,  pour  former  le  développement  de  A2,  il  ne  reste  plus  qu’ 
obtenir  les  développements  des  quantités  r' , r'  cos/',  r'sin/'. 

30.  Posons,  comme  au  n°  19, 

M'  = c'  -+-  n'  t, 

et  soit,  pour  abréger, 


nt  = c -+■  A,  sinco  — As sin 3 a — A, sin 5m  — A, sin  70)—  ... . 


On  aura 

ï\!  Tl f 

M'  = c'  H c (A,  sin — A3sin3co  — A,  sin 5m  — ...  1 . 

n n v 


Posons 

, n' 

c H c — Ci, 

n 


nous  aurons,  tous  calculs  faits, 


M'  =-<?,+  [ i,822o333j  sinw  — [2,836062]  sin 3m  — [3, 181712]  sinÔM 
— [5,4i8]  sin7M  — [7,520]  singM. 


On  a d’ailleurs 


r'  = 1 — ecosw', 
r'cos/'  = cos  — e', 
r'  sin/'  = sji  — e'2  sin 
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3; 


et 


en  négligeant  les  termes  du  quatrième  ordre  par  rapport  à e' . 

Ces  formules,  réduites  en  nombres,  donnent 

r'2=  [o,oooi834]  — [2, 525796] cos M'—  [4,  i486]<?0S2M'—  [6,08]  cos3M', 

r'cosf  — — [2,400873]  + [ 1 ,9999541  ] cosM' 

-i-  [3,92367  | cos2M'+  [4,0237]  cos  3M'-f-  [6,20]  cos  4M', 

r'sin  f'=  [ 1,9999236]  sinM'-l-  [3,923648]  sin  2M' 

-+-  [4,0237]  sin  3M'+  [6,20]  sin  4M'; 

de  sorte  qu’il  ne  reste  plus  qu’à  développer  les  sinus  et  cosinus  de  M' 
et  de  ses  premiers  multiples,  suivant  les  sinus  et  cosinus  des  multiples 
de  co. 


31.  Soit,  pour  abréger 


M'  = c,  -+-  a. 


On  aura 


coscM'=  cos ict  cos ia  — sincc,  sin ioc, 
sin  /M'  = sinic-,  cos  ioc  -f-  coscc,  cos  ioc; 


de  sorte  que,  pour  développer  A2  suivant  les  sinus  et  cosinus  des  mul- 
tiples de  c,-  et  des  multiples  de  w,  il  suffira  de  développer 


On  a trouvé 


a.  — 0,6637939  sin co  — o,o685586  sin  3m 

— 0,0015195  sin5M  — 0,0000262  sin  7M  — o,oooooo3  singe». 
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Posons 


j3  = 0,0015195  sin5co  + 0,0000262  sin7Jo  + o,oooooo3  singco; 
m — 0,663798g, 
n — o,o685586, 
y — m sin  w — « sin  3m. 

Alors 

a = y — [3; 

d’où 

« 3' 

sina  = siny  — (3  cosy  — — siny, 

„ 32 

cos  a = cos  y + (3  siny  — cos  y, 

en  négligeant /33,  ce  qui  est  visiblement  permis. 

Le  calcul  de  /32  n’offre  aucune  difficulté;  on  trouve 

Ÿ (33=  0,0000006  — 0,0000006  cos  ioto. 

Le  développement  de  siny  et  de  cosy  est  plus  difficile;  011  a procédé 
comme  il  suit. 


32.  On  a 


siny  = sin  [ m sin  w]  cos  [n  sin  3m]  — cos[m  sin  w]  sin  [n  sin 3m], 
cos  y = cos  [m  sinM]  cos[«  sin  3m]  -+-  sin  [m  sin  m]  sin  [n  sin  3 m]. 


D’autre  part,  on  a,  par  la  formule  de  Maclaurin, 


. r . n m . m3  . , m3 

sin  m sinM  — — sinM  — r sin3M  „ . - sin'&)  — . . ., 

1 1 1.2.3  1.2. 3. 4. 5 


r . n . , 

cos  msinu  =1 sinJM  H 5 — 7 sin'M  — . . . , 

J 1.2  1 .2.3.4 


développements  qu’il  est  nécessaire  de  pousser  ici  jusqu’à  la  neuvième 
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puissance  de  sinw.  On  a trouvé 


log — — 1,8220333, 
log  — — - = i,343o366, 
log  = 2,6879486, 


log  77^34  = 3’9°79^> 


lot 


m3 


1 .2. 3.4*5 


= 3,o3o985, 


tïl  — 

logrr3T4T5T6  = 4>°’48’' 


log 


/??' 


log 


1 .2. 3. 4-5. 6. 7 
m8 


log 


1 .2.34-5.6.7.8 

ni* 

1 .2.34-5.6.7.8.9 


— 5,o52, 

==  7>97» 

= 8,8- 


On  a employé  les  formules  suivantes  : 

2 sin2«  = — cos  2&j  -I-  1, 

23  sin*  03  = cos4«  — 4 cos  203+3, 

2.isin6o3  = — cos6o3  + 6cos4«  — i5  cos  203  + 10, 

2*  sins03  = cos  8 03  — 8 cos 603  + 28  cos 4 w — 56  cos  2 03  + 35, 

2*sin*o3  = — sin  3 03  + 3 sino3, 

2f  sin3  03  = sin5o3  — 5 sin  3 03  + 10  sin  03, 

26 sin7 03  = — sin  703  + 7 sin  5 03  — 21  sin  3 03  + 35  sin  03, 

28sin9  03  = sin  903  — 9 sin  703  + 36  sin  5 03  — 84  sin  3 03  + 126. 
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L’application  de  ces  formules  a donné 

sin  [m  sino3]  = [ 1,7978896]  sino3 

+ [2,073898]  sin3o3  + [5,819]  sin  5o3  + j 7,3]  sin  703, 

cos  [m  sino3]  = [ 1,9507742] 

+ [1,0259859]  cos 2 03  + [4,9952]  cos4«  + [6,556] cos 6 CO. 
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On  a procédé  de  même  à l’égard  de  sin  [/isin3w]  et  de  cos  [nsin3w]. 
On  a trouvé 

log  — = 2,836o62, 

1 


l°g  ^ = 3,37109, 

lo«,.",.3=5’73°- 


log 


H4 

12.3.4 


7>96> 


et  l’application  des  formules 

2 sin2  3 « = — cos6w  + 1, 

23sin43  w = cos  1 2&>  — 4 cosô co  —h  3, 

22sin33 w = singw  — 3 sin3oi, 

a donné 

sin[rasin3w] ~ [2,8358o66]  sin3w  -+-  [5,1271  J singw, 

cos  [11  sin  3 <so]  = [1,9994895]  +-  [3,06989]  cos6w  -t-  [n\  cosi2w: 

on  en  a conclu 

sin-/  = 0,6235578  sin  w — o,o4g34 1 gsin  3 w — o,oo3g4° 1 sin  5 w -4-  o,ooo3345  sin  7 w 
— o,ooooo5i  singw  — 0,0000007  sin  1 1 w, 
cosy=  0,8921984  + 0,1275555  cos2w-  0,0204610  cos4w  -1-  o, 0006463  cos6w 
-t-  0,0000643  cos8w  — o, ooooo36  cos  iow  + 0,0000002  cos  12 w. 

33.  Des  valeurs  précédentes  on  déduit  par  la  simple  multiplication 

^ sin 7 = — o,ooooo3o  — 0,0000373  cos2w  -+-  0,0004732  cos4<j> 

— o,ooo4656  cos6w  + 0,0000294  cos8w 

-f-  0,0000035  COS  IOW  — 0,0000002  COS  12W, 

j3  cos  y = — 0,0000160  sin  w + 0,0000966  sin  3 w -4-  0,00 1 3574  sin  5 w 

_i_  0,000 1 2o3  sin  7 w — 0,00001 35  singw  -f-  0,0000002  sin  1 1 w, 

ri  sin  y — 0,0000004  sinw  -+-  0,0000002  singw  — 0,0000002  sin  1 1 w, 

Ê1  cosy  = o,oooooo5  -t-  0,0000001  cos2w  — o,oooooo5  cos  iow, 
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et,  en  rassemblant  ces  résultats,  on  trouve 


4< 


siri  a = [ i ,■5948876]  sin  co  — [ 2,694065]  sin  3 6>  — [ 3,72407 1 j sin  5a> 
4-  [4,33o8]  sin  7 &)  4-  [6,914]  singw  — [ 7,8]  sin  1 1 co, 

cos cf.  = j 1 ,g5o45g8]  -4-  [ 1 , io557  17]  cos  2 co  — [2,300765]  cos4&> 

4-  [4,2570]  cos 6 « 4-  [5,972]  cos  8 co  4-  [7,6]  cos  10  co. 


34.  On  a déduit  de  ces  développements  ceux  des  sinus  et  cosinus  des 
multiples  de  u.  par  les  formules 

sin2a=r2sina  cos  a, 
cos2a  = 1 — 2 sin’a, 
sin3cc  = 2sina  cos  2 a -4-  sin  a, 
cos  3 ce  — 2 cosse  cos  2 a — cos  a, 
sin 4«  = 2 sin  2a  cos 2 ce, 
cos4«  — 2 cos*2a  — 1, 

et  l’on  a obtenu  les  résultats  suivants  : 

sin  2 a — [0,01  n 5]  sin 00  4-  [3,4go]  sin 3c.)  — [2,4518]  sin5co 
-4-  [4,87]  sin  7 co  -4-  [4,24]  singe.)  — [5,o  J sin  1 1 co, 

cos2a  = [ 1,78439]  -i-  [ 1 ,653 ig]  cos2co  — [2,7406]  cos4&> 

— [3,6464]  cos6c^  4-  [4,8g]  cos8co  4-  [5,23]  COS  106), 

sin  3 u — [ o,o32  38]  sinco  4-  [1 ,3o68i  ] sin3eo  — [ 2,81 55]  sin5&> 

— [3,387]  sin  7004-  [3,037]  sin 9 o) — [5,3]  sin  1 1 &>, 

cos3a  = [ 1,4021 1]  4-  [ 1,91 100]  cos2co  — [ 2,65g6 j cos4&> 

— [2,3755]  cos6co4-  [3,312]  cos8&) 

-t-  [4,4 1 ] cos  106)  — [5,3]  COS  126), 

sin 4 ^ = [ 1,898]  sin&)  4-  [ 1,708]  sin3&)  — [2,92]  sin5co 

— [2,15]  sin7&)4-  [3,3]  singeo, 

C0s4a  = — [2,72]  4-  [0,0195]  COS  2 6)  4-  [2,81]  cos4&)  — [ 2,78]  C0S6&). 

6 
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35.  Pour  parvenir  à la  formation  des  quantités  rr'  eos/cos/', 
on  a d’abord  obtenu  les  résultats  suivants  : 

r eos/cos  x = — 0,02205 66  4-  0,3359186  cos 2 w 4-  0,0341075  cos4  co 

— 0,0031228  cos6co  — o,oooo4g8cos8o)  o,ooooi75cosioo) 

-t-  o,oooooo5  COS  1 2 co, 

rsin/sin  « = 0,2909391—  o,3o43g56  cos  20)  4- 0,0102/22  cos4co 

-t-  0,0032127  cos6m  -t-  0,0000072  cos8co  — 0,0000057  cos  10 co 

-+-  0,0000001  COS  12  60, 

r cos  / sin  a = — o,  162871 3 sin  co  + o,  1 ig533o  sin3co  — 0,0064  i5g  sin  5 w 

— 0,00 1 2i4g sin 7C0  4-  0,0000129  singco  4-  0,0000020  sin  1 1 co 

— 0,0000002  sin  1 3co, 

r sin/cos  « = 0,7777946  sinoo  4-  0,0994758  sin3co  — o,oo66563  sin5co 

— 0,0002463  sin7co  4-  0,0000410  singco  4-  0,0000018  sin  1 1 co, 

r cos/cos2a  = o,o55i  1 4-  0,20382  cos2co  4-  o,og4g4  cos4«  — 0,00810  cos 6» 

— 0,001 12  cos8co  4-  0,0001 3 cos  ioco, 

/ sin  /sin  2«  = 0,4801 3 — 0,46164  cos 200  — o,o3i8i  cos 4 00  4-  0,01 320  cos6c> 
4-  0,00021  cos8co  — 0,00009  cos  i°«, 

fcos/sin  2 « = — 0,25210  si  no  4-.  0,1 8845  sin  3co  4-  0,00692  sin5co 

— o,oo53g  sin  704-  0,00001  singco  4-  0,00004  sin  1 1 co, 

r sin/cos2«  = 0,36766 sin w 4-  0,2.5728 sin 3co  — o,oi56i  sin5c> 

— o,oo32i  sin  7 co  4-  0,00026  singco  4-  0,00001  sin  z 1 co, 

/•cos/ cos 3 a = o, i445o  4-  o,o5oo6 cos2co  4-  o,i58i4  cos4«  — o,oo338  cos 60 

— 0,00480  C0S8C0  4-  0,00028  COS  IOCO  4-  0,00006  COS  I 2 Ci) , 

/•  sin/ sin3«  = o,5orj5r]  — 0,39197  cos2co  — o,i434gcos4«  4-  o, 02568  cos6co 

4-  0,00268  C0S8&)  — 0,00046  COS  IOCO  — 0,00001  COS  120), 

/ cos / sin  3«  = — 0,227 16  sin  co  + °>  i8o52  sin  3co  4-  0,047 1 5 sin  5 co 

— 0,01 1 3 1 sin 70)  — 0,00082  singco  4-  0,00018  sin  1 1 w 

— 0,00001  sin  i3co, 
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rsin/ cos 3 a = — o, i3o34  sin  6)  H-  o,4i  1 88  sin  3 m 4-  0,00377  sin  5m 

— o,oi258sin7  w -+-  o,ooo4i  singM  4-  o,oooi5sin  1 1 m 

— 0,00001  sin  i3m. 

rcosf  cos4«  = o,2o5  — 0,057  C0S2W  -+-0,186 cos4«  -+-  0,021  cos6m 

— 0,012  cos8m, 

rsin/ sin4«  — 0,379  — o»1 1 9 cos 200  — 0,291  cos4w  -+-  0,023  cos6m 
4-  0,008 cos8m, 

rcos/sin4«  = — 0,099  si n&j  -4-  0,106  sin  3 m -4-  o,  104  sin  5m  — 0,0 12 sin  700 

— o,oo3  singM, 

r sin  / cos 4 « = — 0,521  sine.)  -+-  o,458  sin 3m  4-  0,076 sin  5m  — 0,027  sin  7 m 

— 0,001  singM. 

36.  En  substituant  tous  ces  résultats  dans  la  formule  du  n°  29,  on  a 
finalement  obtenu  le  développement  de  A-  suivant  les  sinus  et  cosinus 
des  multiples  de  c,  et  des  multiples  de  w 

AJ  = U0  -+-  U,  cosc,  -4-  U2  C0S2C,  4-  U3  cos  3 c,  4-  U,  cos  4 c, 

4-  Vi  sine,  4-  V3  sin2c,  4-  V3  sin 3c,  4-  V4  sin 4c, , 

les  quantités  U et  Y ayant  les  valeurs  suivantes  : 

U„  = i,5i54o25  — o,o386io3  sinM  — 0,432.372  cos2m  — o,ooi4o6sin3M 
4-  o,o43223  cos4«  — 0,000028  sin5M  4-  0,002.386  cos6m 

— 0,000001  COS  7 M -I-  0,000074  C0s8m  H-  0,000002  COS  IOM, 

U,  = — o,3i  1223  4-  1,548637  sinM  — 0,048616  C0S2M  — 0,037923  smi3m 

— o, 046991  cos4&)  — 0,00012.3  sin5M  4-  0,000004  cos6m 
4-  0,002002  sin  7 m 4-  o,oooo43  cos 8m  4-  0,000046  singM 

— o,ooooi3  cos  IOM, 

V,  — — 0,439768  — 0,971 178  sinM  — o,o65854  C0S2M  4-  o,o233go  sin  3m 

— 0,074 i3o  cos 4 m — 0,000 1 4 * sin 5m  — 0,00001 5 cos 6m 

— o,ooio4gsin7M  4-  0,000072,  cos8m  — 0,00002.9  singM 

— 0,000020  cos  IOM, 
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U2  = — 0,00481 1 -4-  0,00861 1 si  11  m -4-  0,002148  cos  2co  -4-  0,000878  sin  3o» 

— 0,000674  cos4o)  — o,ooo3i2  sin  5 m — 0,000042.  cos6&) 

-4-  o,oooo32  sin  7 o)  -+-  0,000008  cos8oj  -4-  o,ooooo3  sin 90), 

Vj  — — 0,006825  — 0,007380  si n co  + o,oo3473  C0S2&)  — o,ooo552  sin  3c.) 

— 0,00090 1 cos4o>  -4-  0,000196  sinôc.)  — 0,000067  cos6&) 

— 0,000021  sin  76)-+-  0,0000 1 3 cos 8m  — 0,000002  singe.) , 

U,  = — 0,000073  -4-  0,000018  sine.)  -+-  0,000037  cos2m  + o,oooo3gsin  3o) 

— 0,000002 cos4c<)  — 0,000008 sin5M  — 0,000002  cos6m, 

V.  — — 0,000 1 1 4 — 0,000010  sin  m -4-  o,oooo5g  C0S2M  — 0,00002.5  sin  3 m 

— o,ooooo3  cos4 e»  — o, 000004  cos6m, 

U,  — — 0,000001  — 0,000001  sinM  h-  0,000001  sin3e», 

V4  = — 0,00000 1 . 


37.  Les  premiers  termes  de  la  formule  précédente  sont 


I,5i54o25  — 0,4323720  C0S2M, 

— sine,  [0,439768  -4-  0,971 1 78 sin c.j] , 

— cos c,  [o,3i  1223  — 1,5486368  sin  m]. 

A la  simple  inspection  de  ces  termes,  on  voit  que  la  méthode  de  M.  Gyl- 
den  11e  peut  être  appliquée  : il  faut  diviser  l’orbite  en  plus  de  deux 
parties. 

On  peut  placer  un  point  de  division  au  périhélie  en  posant 


sinM  = 


sm2 — 
2 


I — COSC.), 
2 


5 


sin  w — — sin2 


0)  2 

2 


I — COS  «2 
2 


En  adoptant  la  première  substitution  et  faisant  varier  m,  de  zéro  à n,  011 
aura  la  partie  de  l’orbite  comprise  entre  le  périhélie  et  le  point  (a, , r,). 
La  seconde  substitution  se  rapporte  à la  partie  symétrique  par  rapport 
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au  grand  axe.  Appliquant  la  première  substitution,  on  trouve 

A2  = 1,4072095  — 0,925357  sine,  -+-  0,463095  cosc, 

— 0,432372  cosw,  + 0,485589  cosw,  sin  c , — 0,7743184  cosw,  cosc, 

•4- 


La  possibilité  d’appliquer  la  méthode  de  M.  Gylden  commence  à appa- 
raître; en  toute  rigueur,  on  pourrait  le  faire;  les  termes  indépendants 
de  a),  sont  les  plus  grands.  Néanmoins,  il  paraît  convenable  de  diviser 
encore  cette  partie  de  l’orbite.  On  peut  le  faire  en  posant,  pour  la  partie 
la  plus  voisine  du  périhélie, 


COSW, 


cos2  ~ w', 


I -h  cosw', 
2 


et,  pour  l’autre  partie, 

cosw,  = — COS2  4 w", 


1 -t-  cosw" 
2 


w"f  variant  de  — n à zéro  et  w',  de  zéro  à tt. 

Si  l’on  adopte  la  dernière  formule,  les  ternies  donnés  précédemment 
deviennent 

1,623395  — 1 , 1 68 1 5 1 sine,  + o,85o254  cosc, 

-1-  0,216186  cosw"  — 0,242794  sine,  cosw"  -+-  0,387  '^9 cosc,  cosw" . 

L’argument  w"  se  rapporte  à la  partie  de  l’orbite  la  plus  voisine  de 
l’orbite  de  la  Terre,  et,  par  suite,  à celle  où  les  perturbations  seront 
les  plus  grandes.  Il  est  visible  que  toute  la  difficulté  de  développe- 
ment des  puissances  négatives  de  A résidera  dans  celui  des  puissances 
du  polynôme. 

i,6233g5  — 1,1 68 1 5 1 sine.  H-  o,85o254  cosc,. 

L’introduction  de  l’argument  elliptique,  d’après  les  indications  don- 
nées au  n°  20,  rendra  ce  développement  très-facile. 


38.  Hansen,  dans  le  Mémoire  cité  sur  le  calcul  des  perturbations  des 
comètes,  dit  que  le  développement  général  de  A-3,  dans  le  cas  parti- 
culier que  nous  avons  traité,  serait  impraticable.  Aussi  a-t-il  dû  donner 
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à la  variable  « des  valeurs  particulières,  développer  les  puissances  de 
.A-3  pour  toutes  ces  valeurs  particulières  suivant  les  sinus  et  cosinus  des 
multiples  de  c{  et,  par  des  quadratures,  reconstituer,  au  moyen  de  ces 
développements  particuliers,  le  développement  général  procédant  à la 
l’ois  suivant  les  sinus  et  les  cosinus  des  multiples  des  deux  arguments. 
Comme  on  l’a  vu  au  dernier  paragraphe,  la  méthode  de  M.  Gylden  ne 
s’applique  pas  'a  la  portion  d’orbite  que  Hansen  avait  embrassée  dans 
son  calcul  ; il  faut  la  subdiviser  en  plusieurs  parties.  Malgré  la  longueur 
des  calculs  résultant  de  cette  subdivision,  il  semble  que  les  idées  de 
M.  Gylden  aient  un  avantage  considérable.  La  marche  à suivre  est  à la 
fois  très-claire  et  très-simple;  les  diverses  portions  de  l’orbite  peuvent 
être  calculées  séparément  par  plusieurs  personnes,  tandis  que  la  mé- 
thode de  quadrature  employée  par  Hansen,  méthode  extrêmement 
longue  d’ailleurs,  n’a  aucun  de  ces  avantages. 


Vu  et  approuvé. 


Paris,  le  24  mai  1876. 

Le  Doyen  de  la  Faculté  des  Sciences, 


MILNE  EDWARDS. 


Permis  d'imprimer. 

Le  Vice-Recteur  de  l’Académie  de  Paris, 
A.  MOURIER. 


SECONDE  THESE. 


PROPOSITIONS  DONNÉES  PAR  LA  FACULTÉ. 


Résolution  des  équations  du  cinquième  et  du  sixième  degré  au 
moyen  des  fonctions  elliptiques. 


Vu  et  approuvé. 

Paris,  le  24  niai  1876. 

Le  Doyen  de  la  Faculté  des  Sciences, 

MILNE  EDWARDS. 


Permis  d’imprimer. 


Le  Vice-Recteur  de  l’Académie  de  Paris, 


A.  MOURIER. 


2<J4ü  paris 
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